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Pauli Equation on a Curved Surface and Rashba Splitting on a Corrugated Surface 
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The Schrodinger equation for a spinless charged particle on a curved surface under an elec- 
tromagnetic field has been obtained by adopting a proper gauge which allows the separation 
of the on-surface and transverse dynamics. [Phys. Rev. Lett. 100 (2008) 230403] As its exten- 
sion, I provide the Pauli equation for a charged spin- 1/2 particle confined to a curved surface 
under an electromagnetic field. Energy spectra of a sphere and a corrugated surface to which a 
particle is confined are given as simple applications of the equation. The energy levels obtained 
exhibit splittings due to the relativistic effect known as the Rashba effect. 
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Much attention has been paid to low-dimensional sys- 
tems because they show exotic phenomena, 1-8 essentially 
distinct from solid state materials. Recent development 
of nanotechnology facilitates construction of nanostruc- 
tures with curved geometry 9-14 and those systems can 
be experimental platforms for study of the interesting 
phenomena. Curved geometry induces the geometric po- 
tential, which affects the dynamics of an electron moving 
on the curved surface, even when an electrostatic poten- 
tial is absent. 

The Schrodinger equations on curved surfaces have 
been used as tools for theoretical investigations on such 
systems. Their formulations have been based on the two 
methods. The one was proposed by DeWitt 15 and the 
other by da Costa. 16 While the former method regards 
a curved surface to be a fully two-dimensional space and 
starts from a Lagrangian for a curved space, the latter 
method regards the curved surface as a two-dimensional 
system embedded in the flat three-dimensional space. 
The method established by da Costa, which is called 
the thin-layer method, has been widely used despite its 
impossibility of inclusion of arbitrarily oriented mag- 
netic fields. Da Costa also provided the one-dimensional 
Schrodinger equation for a curved thin tube. Along 
this line, Takagi et al. 17 studied the geometry-induced 
Aharonov-Bohm effect by taking into account the torsion 
of a tube. 

Ferrari et al. 18 recently adopted the thin-layer ap- 
proach and rigorously demonstrated by choosing a 
proper gauge that the separation of the on-surface and 
transverse dynamics under an electromagnetic field is 
possible without approximations. Theoretical investiga- 
tions on nanostructures are thus expected to treat more 
various situations in future. 

It should be pointed out here that there also exists 
another rich field of studies for formulation of quantum 
mechanics on a more general constrained system by em- 
ploying the possibility of incquivalcnt quantizations (e.g. 
Refs. 19-22). Ohnuki et al. 21 analyzed quantum mechan- 
ics on S D embedded in (D + l)-dimensional flat space by 
setting up appropriate fundamental algebra, with gauge 
potentials induced. 



Parallel to the extension of the possibility of geomet- 
ric arrangement, the role of relativistic effects such as 
spin-orbit interaction has been increasingly evoking in- 
terests in material physics due to the recent rise of spin- 
tronics. The relativistic effects and the spin degree of 
freedom for an electronic system naturally emerge in the 
Dirac equation. 23 Since the upper two components of a 
four-component spinor are much larger than the lower 
two components in ordinary condensed matter physics, 
the large part of relativistic electronic structure calcu- 
lations have been performed using two-component wave 
functions. The most tractable tool for a quantum me- 
chanical analysis on an electronic system including mag- 
netic properties and relativistic effects is hence the Pauli 
equation, which is for a charged spin-1/2 particle with a 
nonzero mass and contains the lowest-order relativistic 
correction term. It will provide deep insights into rich 
physics on nanosystcms. 

In this Letter I provide the Pauli equation for a particle 
confined to a curved surface under an electromagnetic 
field as an extension of the work done by Ferrari et al. 18 
I reconcile ourselves to two assumptions mentioned below 
for performing variable separation in the present study. 
Excepting those assumptions, the derivation of the Pauli 
equation will proceed on the same strictness as for the 
Schrodinger case. 

The expansion of the Dirac equation for a spin-1/2 
particle of its mass m and charge Q under an electromag- 
netic field using the Foldy-Wouthuyscn method 24 leads 
to the Pauli equation = Hpyj for the upper two com- 
ponents tjj of the four-component spinor. Hp is the Pauli 
Hamiltonian, which neglects the mass term mc 2 and the 
terms on the orders higher than to -2 , given by 

V-®-SB 



4m 2 c 2 



■SxE-U]. 



(1) 



II = — iS/ — — A is the canonical momentum operator and 
the magnetic field B = V x A is the rotation of the vector 
potential. S = cr/2 is the spin operator and cr is the 
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Pauli matrix. The contribution from the divergence of the 
electric field V • E, called the Darwin term, is absorbed 
into the scalar potential V. All spin-orbit interactions 
which will appear below come from the last term on the 
right hand side of eq. (1). Let us consider a coordinate 
transformation from the Cartesian coordinates xi a ) {a = 
x, y, z) to the curvilinear coordinates qi(i = 1,2,3). The 
letter in a parenthesis stands for the Cartesian coordinate 
and not in a parenthesis for the new coordinate. The 
dreibein field and its inverse are defined as 

(2) 



(a) 9x {a 



which satisfy the conditions e\ e 



-(a) ~ "> ' c (b) ~ °(b) 

due to the chain rule of derivative. Summation is im- 
plied over the repeated index. The metric tensor in the 



new coordinate system is given by Gi 



We 



define v; 



Ha) 



for a vector v 



W(o)e(«) 



where e 



(o) 



is the unit vector along a direction in the Cartesian co- 
ordinate system, from which it follows that V( a ) = e^Wj. 
Using the Laplacian for the curvilinear coordinate sys- 
tem 25 and the covariant derivative D{ = di — the 
Pauli Hamiltonian, eq. (1), is rewritten as 



(o)-»(a) 



H ^=-Li& DiVUGijDj+v -^' 



(3) 



where G = detG,-., and G 1 ^ is the inverse metric tensor. 



We have defined B = B 



i dB 

Amc dt 



using the Maxwell's 



equation V x E = —-^M-- We have defined the matrix 



-h? 1 , where ei abc \ is the Levi- 



tt -> V = V 



fe« = £(abc)0{a)e\ b) e\ c) 
Civita symbol. The Pauli equation is obviously invariant 
under the following gauge transformation with an arbi- 
trary scalar function 7: 18 

c dt' 
At -> A\ = Ai + dij, 

ip -> ip' = e iQl/c ip. (4) 

Writing explicitly the spatial derivative with respect to 
the new coordinates, we write down the Pauli equation 

as 



dip 
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2m 



1 



Vg 



diiVGG^djip) 
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unit normal vector of S at r 5. It obviously follows that 



-—G ij Aidiib - %G^A l A 1 ip 
c c l 



+Vip- 
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From here the curved surface S confining the particle 
is considered. Let us adopt a coordinate transforma- 
tion such that S is described as 93 = and an arbi- 
trary point r immediately close to the point rj on S 
is given by r(q 1 ,q 2l q 3 ) = T\s(gi,<? 2 ) +?3w(gi,g 2 ), where 
e a = ^f-( a = lj2) and n = e\ x e 2 /|ei x e 2 | is the 



» 



l,2),e 3 a) 



n ■ ei a \ on S. The two- 



dimensional induced metric tensor g a b = e a • e& is con- 
nected with the three-dimensional one via the following 
relation: 

Gab = 9ab + [ctg + t (ag)}abq3 + {ag t a) ab ql, 

G a3 = G 3a = 0, G 33 = 1. (6) 

a a b is the Weingarten matrix, 16 which satisfies ^p- = 
a ab e b . As is done in the case of the Schrbdingcr equa- 
tion, 5, 16, 18 we put the wave function in the form 



x(qi,q2, q 3 ,t) 



(7) 



\J\ + Tr aq 3 + det aq 2 

which ensures the norm conservation condition 
J\ip\ 2 VGd 3 q = J\x\ 2 yfgd 3 q. Substituting eq. (7) 
into eq. (5) and taking the limit 93 — > 0, we obtain the 
equation for X- 



.dx 
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2m 



da(Vgg ab dtx) 



V9 



d a (^g ab A b ) X 



2 -l9-g^A a d bX -%{g ab A a A b + Al) x 
c cr 



iQ, 
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2iQ 
c 



Vsx 
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2mc 



V(a)B( a )X 



Am 2 c 2 



h^E i D j --TraK ,a E i \ X , 



(8) 



where Vs (91,(72) = — -^{(^ a Y / ^ ~ det a] is the well- 
known geometric potential. The electromagnetic field 
is evaluated at g 3 = 0. We here apply a gauge 
transformation with a scalar function 7(91,(72,(73) = 
— J 93 A 3 {qi, 92, z)dz, as done in the Schrodinger case, 18 
so that A' 3 and d 3 A' 3 vanish and A\ and A2 are un- 
changed on S. We divide the scalar potential into the 
on-surface electric part and the confinement part as 
V = V c \ (qi, 92) + V\(q 3 ). The gauge-transformed equa- 
tion reads, from eq. (8), 



■ dx 
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2m 



±d a (^g ab d bX ) - —-^d a (Jgg ab A b ) X 
V9 c V9 



^g ab A a d b x-%g ab A a A bX 
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+ V S X 
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where 



2mc 



<7 (a)-B(a)X + 



iQ 
4m 2 c : 



■hPEiDM + Hsx, (9) 



Ha 



1 



hTEiida-Traffl + Vx (10) 



2m J ' 4m 2 c 2 
is the Hamiltonian describing the dynamics along 93 di- 
rection. H 3 contains 91, 92 and t as parameters. The term 
involving the first-order derivative on the right hand side 
couples the on-surface and transverse dynamics as a rcl- 
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ativistic effect. In the case of the Schrodinger equation, 
this term were absent and hence the variable separation 
would have been done. Let the solution of an eigenvalue 
problem H 3 (f> = £</> of the form 



<f>(Qi,Q2,Q3,t) = f(qi,q2,q3,t) 



(11) 



The condition for a nontrivial solution is then, from eq 
(10), 



QE 

4m 2 c 2 



Tr a 



1/ 







(12) 



with g-f and g± arbitrary, where E( a 
and E = 



We try the form 

f(Qi,Q2,q3,t) = ex P( ± j^Eq 3 )f(q 3 ) and substitute it 
into eq. (12). The condition then becomes 



E tv) 



(z) 



2m 3 A 32?ti3c 4 + 8m 2 c 2 ra 



/ = £/■ (13) 



It implies that the eigenvalue is expressed as 

. Q 2 E 2 iQE 

T ~ — — Tra, 



32m 3 c 4 8m 2 c 2 



(14) 



where e is an eigenvalue for E = 0, independent of qi, q2 
and t. Wc neglect the second term on the right hand side 
of eq. (14) since it is on the order of m -3 . In addition, wc 
take the average of the positive and negative signatures 
of the third term so that it vanishes. The first assumption 
needed for the present derivation is the validity of this 
treatment. The confinement potential is so strong that 
the wave function at any point on S is, as the second 
assumption, expected to be written as 



x{Qi,Q2,Q3,t) = exp ± 



iQE 
Amc 2 



q3\fo(q3)xs(qi,q2,t), 

(15) 



where /o is the nondegenerate solution of eq. (13) with 
the lowest eigenvalue Eq. This assumption is feasible 
as long as low-energy excitation is discussed. The two 
assumptions introduced above are not needed for the 
Schrodinger case 18 thanks to the absence of the rela- 
tivistic term, for which the differential equation of the 
transverse dynamics remains unsolved after the variable 
separation. Hereafter we set eq = 0. Substitution of eq. 
(15) into eq. (9) leads to the equation to be satisfied by 
the surface wave function xs : 



■ dxt 



Ot 



1 

2m 



— d a (^g ab d bXS ) 



c c 



-9lg«b AaAbXs 



+ (Vs + V el ) X s + H sp - relX s, (16) 



where 



H, 



sp— rel 
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■ — 

2mc 
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Am 2 c 



thTEiDa. (17) 



It has been demonstrated that the variable separation is 
possible with the two assumptions for the Pauli equation 
on a curved surface under an electromagnetic field. The 
Hamiltonian of the resultant equation, eq. (16), is simply 
that for the Schrodinger equation on a curved surface 
obtained by Ferrari et al. 18 plus H sp - Te \. The new term 
acts on the two components of the surface wave function 
differently in general and hence physics which is absent 
in the Schrodinger case can emerge. 

In what follows, energy spectra of a charged spin- 
1/2 particle on a sphere and a corrugated surface 
are provided as instructive applications of the surface 
Pauli equation. Analyses of their spinless and nonrela- 
tivistic electronic properties have been done using the 
Schrodinger equaiton for curved surfaces. 1 ' 2 ' 6 

As the first example, let us calculate the energy spec- 
trum of a charged spin- 1/2 particle confined to a sphere 
of radius R. We here adopt the spherical coordinates as 
a new coordinate system. We assume that there is no 
magnetic field and a static electric field with a constant 
amplitude penetrates perpendicularly to the sphere: E = 



Ee r . In this case Vs 
and 



Vol 



0, E r — E, Eg 



E,i 







COS( 



h r4> 



cot 9 cos ( 
R 

cosci 



R 



T( X ) 



(x) 



R 

cot 6 sin ( 
R 



-a 



(y) 

i> 

- a (y) 

cote 



(18) 



on the sphere. The Pauli Hamiltonian on the sphere is 
thus 



H 



1 



2mi? 2 



■ — -do sin Ode H o 



sin ( 



ia ( d<f, —ie ^(8$ — icotBd^) 

R [ie^idg +icot0d tj> ) -d+ 

L 2 



2a 



LS, 



(19) 



2mi? 2 R 

where a = — ^J?^ and L is the orbital angular momen- 
tum operator. This Hamiltonian contains the spin-orbit 
interaction in the form well known in condensed matter 
physics, which is absent in the Schrodinger case. Spinor 
spherical harmonics 26 




(20) 



is a simultaneous eigenfunction of J 2 , J z , L 2 and S 2 with 
the eigenvalues j(j + 1), j z , 1(1 + 1) and 3/4, respectively. 
J = L+S is the total angular momentum operator and j 
can take only / ±1/2. It is easily confirmed that /Ris 
a normalized eigenfunction of the Hamiltonian, eq. (19), 
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with the eigenvalue 
1(1 + 1) 



en 



+ -[j(j + l)-^ + l)-3/4]. (21) 



" J 2mR 2 ' R l 
The energy gap between two levels with common I and 
different j is -^(2/4-1), which is a relativistic effect. Since 
j ± 1/2 = I « mvii = /ci?, where w is the velocity of the 
particle and k is its wave number, in the limit of R — > oo 

k 2 

the energy eigenvalue converges to Eji ~ #- ± ak. This 
energy dispersion is the same as that for the Rashba 
Hamiltonian 27 Hn — 4 — h ctcr • e z x fc, which describes 
a charged particle moving freely on xy plane under an 
electric field along z axis, leading to the spin splitting as 
a relativistic effect due to the lack of inversion symmetry. 

As the second example, let us calculate the energy 
spectrum of a charged spin- 1/2 particle confined to a 
currugated surface. We assume here that the corru- 
gation is along x direction represented by z = f(x) 
and there is no magnetic field and a static electric 
field with a constant amplitude penetrates perpendic- 
ularly to the surface: E = E(-f'e x + e z )/y/\ + f' 2 
and V e \ = 0. Wc define the new coordinates on S as 
qi = j^dx'i/i + f'{x') 2 ,q 2 = y,q 3 = 0. qi is the line 
length along S. e 1 = (e x + f'e z )/y/l + f 2 , e 2 = e y 

and thus g ab = S ab ,n = (-f'e x + e z )/y/l + f' 2 and 
i f" 2 

— ~ §m (i+f' 2 )' i • nonzcro components of the 



dreibein on S are eL^ 



=3 _ 
-(*) 



"(x) 



JL 



: (y) 



l,ej 



0) 



and 



Vi+7 72 



Using E l = E 2 = 0,E 3 = E,V cl = 



h 12 = 



(.'■) 



K 2i = 



a (x) + f l(J {z) 



h 31 =a 



(y) 



(22) 



on S, we obtain the following Pauli Hamiltonian for the 
surface wave function: 

1 



H=- — {dl + d 2 2 ) + V s 

. ( p. o-(x) + f l(J (z) a 
-ia \ cr (y) di „,„ o 2 



12 



(23) 



where the term proportional to a is responsible for rela- 
tivistic effects. We examine periodic corrugation on S in 
what follows, putting f(x) = hg(x), where the constant 
h measures the height of the corrugation and g is a pe- 
riodic function. When the corrugation is absent (h = 0), 
the system reduces to the ordinary Rashba system 27 and 
the time-independent Pauli equation Hu = eu has two 
plane wave solutions: 

ip, e i{kiqi+k 2 g2) 



2m 



± ak, K 



V2 



1 

±ie l 



(24) 



where k\ and k 2 are the wave numbers, k = \/k\ + k 2 
and Ok is the argument of fci + ik 2 on the complex 
plane. From the perturbation theory, the perturbed en- 



ergy spectrum to the lowest order for small corrugation 
is evaluated as the expectation value of the perturbed 
Hamiltonian with respect to the unperturbed wave func- 
tions. Expansion of the Pauli Hamiltonian, eq. (23), in 
terms of h leads to the lowest-order perturbation 

H {1) = iahg'a (z) d 2 (25) 

and the second lowest-order perturbation 

(2) h 2 g" 2 .ah 2 g' 2 

Since (uj_|i?^ 1 - ) |uj.) = due to the periodicity of g and 



z \H^\ul 



h 2 (g" 2 ) 



T f h 2 (g' 2 )k sin 2 k , the per- 



turbed energy dispersion is given by 



e± = — ± ak I 1 
2m I 



2 / 8m 



which indicates that any small periodic corrugation de- 
creases the Rashba splitting. 

In conclusion, I have provided the Pauli equation for 
a charged spin-1/2 particle confined to a curved surface 
under an electromagnetic field by performing the vari- 
able separation on the manner similar to that for the 
Schrodingcr equation. Energy spectra of a sphere and a 
corrugated surface to which a particle is confined were 
calculated. It was demonstrated that any small periodic 
corrugation decreases the Rashba splitting on a surface. 
The basic equation obtained will be a powerful tool for 
analyses of physics on curved surfaces and nanostruc- 
tures where spin degree of freedom and / or spin-orbit in- 
teractions play important roles. 
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